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Abstract 

Using hyperbolic form convolution with doubly isometry-invariant kernels, the explicit ex- 
pression of the inverse of the de Rham laplacian A acting on m-forms in the Poincare space H™ 
is found. Also, by means of some estimates for hyperbolic singular integrals, L p -estimates for 
the Riesz transforms V*A _1 ,z < 2, in a range of p depending on m,n are obtained. Finally, 
using these, it is shown that A defines topological isomorphisms in a scale of Sobolev spaces 
H* np {W l ) in case m^,f. 
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1 Statement of results and preliminaries 

1.1. The main object of study in this paper is the Hodge-de Rham Laplacian A acting on 
m-forms in the Poincare hyperbolic space (M n ,g). The aim is to prove that A defines topo- 
logical isomorphisms in a range -ff^j )P (lrII n ) of Sobolev spaces of forms defined as follows. For 
< m < n, 1 < p < oo and s £ N, the Sobolev space H^^W 1 ) is the completion of the space 
T> m (W a ) of smooth m-forms with compact support with respect the norm 



\\v\\p,s = J^W^vWp- 

i=0 

Here V® means the i-th covariant differential of r], and for a covariant tensor a 

\a(x)\ p d/j,(x) 



\a\ being the pointwise norm of a with respect the metric g and d[i the volume-invariant measure 
on HP given by g. The space (M n ) can be alternatively defined in terms of weak derivatives. 
The main result of this paper is: 
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Theorem A. A is a topological isomorphism from H^ p (M n ) to for p 6 (pi,P2) with 

2(n - 1) 11, 

n — 2 + |n — 2m| pi P2 

in case m ^ 2. 

In the exceptional case m = A is one to one but is not a topological isomorphism for 

any p. For this case we obtain as well some weighted estimates. If m = t|, A is known to have a 
non-trivial kernel. Of course, Sobolev spaces H^ np can be considered for non integer s as well, and 
the same results holds by interpolation. 

Notice that the Hodge star operator * establishes an isometry from H^^iW 1 ) to H^_ m ^ p (M n ) 
which commutes with A, and this is why the range {pi ) p 2 ) depends only on \n — 2m\. Notice too 
that the range (pi,P2) always contains p = 2 in the non-critical case |n — 2m| > 1 and that for 
functions (m = 0), the range of p is (l,oo) (see comments below). We point out that the range 
(pi)P2) equals || — || < ^j-, where \i denotes the greatest lower bound for the spectrum of A in 
H^ 2 (M n ), whose value (P) is ft = (n ~ 1 ~ 2m)2 (for m < §). 

For the Sobolev spaces for p = 2, H^ l2 QH n )j another proof of the theorem, based in energy 
methods and valid for an arbitrary complete hyperbolic manifold, is given in BG . The motivation 
for the theorem, as with |BG| . comes from mathematical physics, where most operators exhibit 
A as their principal part and results like the above become essential to establish existence and 
uniqueness theorems. 

Our method of proof is simply to construct an explicit inverse L for A on T> m (M n ) and show 
that there is a gain of two covariant derivatives 

\\Lri\\p tS+2 < const \\rj\\ P;S . 

Thus Lrj plays the role of the classical Riesz transform in the Euclidean setting. The most delicate 
part is of course 

||V (2) Lt7|| p < const \\rj\\ p , Pl < p < p 2 , rj £ V m (M n ) . 

Riesz-type operators such as VA~5 , V^A -1 have extensively been studied in different contexts, 
for the case of functions. On symmetric spaces, they are bounded in L p , 1 < p < oo and of weak 
type (1, 1). This was shown in |AL| for the first order ones in some spaces, and later extended to 
all symmetric spaces in 0. The L p -boundedness holds as well for higher order Riesz transforms 
in symmetric spaces, but not generally the weak type (1, 1) estimate. In more general contexts, 
this has been shown in |Llj . |L2j . |L3j . among others. In case of m-forms, < m < n, as far as we 
know, there are much less known results, and is for those that our result is new. In |Plj . |P2j some 
aspects of harmonic analysis of forms are developed; in particular, the exact expression for the 
heat kernel is given, and it is very likely that from it one can get as well an explicit expression for 
A -1 . Strictly speaking, to prove the result an exact expression of A -1 is not needed, it is enough 
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having estimates for the resolvent both local and at infinity. In |L3| . estimates of this kind are 
obtained and applied to Sobolev-type inequalities for forms, and they might work for this purpose 
too. However, we feel that our approach, that we next describe, is more elementary and might be 
interesting in itself. 

The de Rham Laplacian A is invariant by all isometries 92 of EP. These form a a group that 
we denote here by Iso(HP). Denoting by ip*{rj){x) = r](<p(x)) the pull-back of a form 7] by <p, this 
means that A and ip* commute, for all ip G Iso(H n ). Therefore the inverse L of A should commute 
too with Iso(HP). Among all isometries of HP, the hyperbolic translations Tr(EI n ) constitute a 
(noncommutative) subgroup, in one to one correspondence with H" itself. In section 2 we do some 
harmonic analysis for forms in HP and introduce hyperbolic convolution of forms to describe all 
operators acting on m-forms and commuting with Tr(EI n ). In a second step (subsection 2.2) we 
characterize the hyperbolic convolution kernels k(x,y) corresponding to operators commuting with 
the full group Iso(H n ). 

Once the general expression of an operator commuting with Iso(H ra ) has been found, we look 
for our L among these. This corresponds to L having a kernel k(x,y) which is a fundamental 
solution of A in a certain sense, and having the best decay of infinity. This kernel turns out to be 
unique for m 7^ %, we call it the Riesz kernel for m-forms in HP, it is found in subsection 3.1 
and estimated in subsection 3.2. Section 4 is devoted to the proof of the ZP-estimates. Here we 
use standard techniques in real analysis (Haussdorf- Young inequalities, Schur's lemma, etc.). For 
the second-order Riesz transform, to show its boundedness in the specified range (pi,P2) needs 
considering some notion of "hyperbolic singular integral" . There exist some references dealing with 
this, e.g. |L4j . [I], and giving some criteria for L p -boundedness that might apply; however, as 
the singular integral arises locally we have found easier and more elementary to treat it with the 
classical Euclidean C alder on- Zygmund theory as a local model, and patch it in a suitable way to 
infinity. 



1.2. We collect here several notations and known facts about EP. We will use both the unit ball 
model B n with metric g = 4(1 — | a? | 2 ) 2 ^ i dx % dx % and the half-space model W\ = {x n > 0} with 
metric g = x~ 2 ^ dx l dx % . Both models are connected via the Cayley transform tp : — > M n given 
in coordinates by 



E 



Vi = — 1 , i = l,..., n-1; Vn = — l • 

5>2 + (x„ + l) 2 J>2 + (x n + l) 2 

i 1 

We denote by e G HP the point (0, 0, . . . , 1) G Wl or G B n . 

The metric g defines a pointwise inner product (a,/3)(x) between forms at x, for every x G HP, 
and a volume measure d/i. In the ball model dfi is written 
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d/j,(x) = 2 n (l — |x| 2 ) n dx 1 . . . dx 11 , and dfi(x) = x^dx 1 . . . dx n in the half-space model. We denote 
by ( , ) the pairing between forms that makes 2 (HP) a Hilbert space 



{a, (3) = / (a, (3)(x) dn(x) . 



We write \a\ and ||a|| for the pointwise and global norms, respectively, of the form a. In terms of 
the Hodge star operator * the inner product can be written too 



(a, (3) = f 



a A * (3 . 



The group Tr(H n ) of hyperbolic translations is in one to one correspondence x i— > T x with H" 
through the equation T x (e) = x. The equations of z = T x y are better described in the half-space 
model by 

Z{ — X n yi + Xi , i — 1, ... ,71 1, z n — x n y n . 

It is easily checked that indeed Tr(H n ) is a (non-commutative) group. The inverse transformation 
of T x will be denoted S x . Another explicit isometry ip x mapping e to x, satisfying ip~ 1 = ip x is 
given in the ball model by 

, x (\x\ 2 - l)y + (\y\ 2 - 2xy + l)x 
MV) = \xW-2xy+l • (1) 

Since the isotropy group of is the orthogonal group O(n), the general expresion of ip € Iso(EI ra ) is 
V = tPx ° U, with x = ip(0). 

The hyperbolic (or geodesic) distance between x,y £ M n is written d(x,y). We will rather use 
the pseudohyperbolic distance r = r(x,y), related to d by the formula d(x,y) = 2arctanhr(a;, y). 
The explicit expression of r(x,y) 2 in the M" model and the B n model is respectively 

r2= i tri ■ 

\x - y\ z + 4x n y n 

r 2 = \cp x (y)\ 2 = \ X ~ x,y £M n . (2) 

(1 — |x| z )(l — \y\ ) + \x — y\ z 

Associated to the group of translations we have the basis of orthonormal translation-invariant 
vector fields Xi(x) = (T x )*(Xi(e)), such that Aj(e) = They satisfy Xi(u o T x ) = (X.iu) o T x 
for every smooth function u. We will denote by w l {x) the dual basis of Xj, which accordingly is 
orthonormal and tranlation invariant too: T*w l = w l . Their expression in the R™ model is simply 

Xi(x) = x n -^- , w l (x)=x~ l dx\ i = l,...,n. 

Because of their translation-invariance property, the (Xi,w l ) are more suitable than the (Xi,rf) 
defined in the ball model B n by 

Yi{x) = (1 ~ 2 |X|2) ^ , = 2(1 - Ixl 2 )" 1 ^ . 
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For an increasing multindex / of length |/| = m we write w 1 = w H A w 12 A • • • A w %m , and similarly 
dx 1 or rj 1 . The {w J }j is an orthonormal translation-invariant basis of m-forms. 

Recall that the de Rham Laplacian is defined A = dS + 5d, where 5 is the adjoint of d with 
respecto to ( , ). Although strictly speaking not needed, the following expression of A in w 1 - 
coordinates will simplify the analysis at some points. If a = J2i a i w ^ a computation shows that 
in case n J 

(Aa)j = Aaj + 2 X k a Jk - p(n - p - l)aj . (3.1) 
Here Jk means the multindex obtained replacing k by n. In case n £ J, 

(Aa)j = Aaj-2^2 Xwu - (1 - p)(p - n)aj (3.2) 

where IJ means the multiindex obtained replacing n by I. For a function / 

n 

Af = -Y,Xff + (n-l)X n f. 
i=l 

In the ball model, with usual coordinates, 

A / = -I(l-|x| 2 ) 2 + (3 .3) 

2 Translation invariant and isometry invariant operators on forms 

2.1. We are interested in finding the general expression of an operator acting on m-forms, and 
isometry- invariant. In a first step we consider translation-invariant operators acting an m-forms; 
these are described by what we might call hyperbolic convolution as follows. Let k(x, y) be a double 
m-form in x, y and define 

(C k a)(x)= I a(y) A* y k(x,y) = (a,k(x,-)) , a G V m {W l ) . 

If T z is a translation with inverse S z 

C k {T* z a){x) = [ (T* z a)(y) A * y k(x,y) = [ a{T z y) A * y k(x,y) 




a(y) A * y k(x,S z y) 



T* z {C k a){x) = C k a(T z x) = I a(y) A * y k(T z x,y) . 
Therefore C k is translation invariant if k is doubly translation invariant in the sense that 

k(x,y) = k(S z x,S z y) VS Z . 
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Using the translation-invariant basis of m-forms w 1 we see that the general expression of k is 

H x > y) = ^2 ki >j( x > yW( x ) ® wJ {y) 
i, j 

with kj(x, y) doubly-invariant functions, that is, of the form kjj(x, y) = ai t j(S y x) for some function 
(or distribution) ai t j. If Sq denotes the Delta-mass at e and 

S(x, y)=Y^ S (S y x)w I {x) w J {y) 
I, J 

then formally 

a(x) = / a(y) A * y 5(x,y) . 



If P is an operator on m- forms commuting with the T y ,S y , we will thus have 

Pa(x)= a(y) A* y P x (5(x,y)) 



and indeed k(x,y) = P x (5(x,y)) is formally doubly- invariant. This shows, in loose terms, that 
the operator Ck of convolution with a doubly translation invariant kernel k gives the general 
translation-invariant operator acting on m-forms. If 



k(x, y) = ^2 ai tJ {Syx)w 1 ' (x) w J (y) 



i.J 



and a(x) = aj(x)w I (x), then C^a has in the basis w 1 (x) coefficients given by 

(C k a) I (x) = Y] / ai i j(S y x)aj(y) dfj,(y) . 
j Jw n 

Thus in the basis w 1 everything reduces of course to convolution of functions. For a function 
convolution kernel a{S y x) and a test function u £ T>(M n ) we may think in 



C a u(x) = / u(y)a(S y x) dfi(y) 

as an infinite linear combination of inverse translates a(S y x) of a(x). Since the vector fields Xi 
commute with translations, it follows that whenever everything makes sense 

X i {C a u)=C Xl aU. (4) 

We point out that this convolution is not commutative, C a u is in general different from C u a. 
Correspondingly, XiC a u — C a XiU is in general not zero; in fact one can easily show ( BG ( lemma 
3.1]) that these commutators are linear combinations of other convolution operators built from 
a(Syx). 
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2.2. Let P be a generic translation-invariant operator acting on m-forms. We have seen in the pre- 
vious subsection that we can associate to P a doubly-translation invariant kernel k(x, y) so that P = 
Cfc. By the same argument as before, P will be isometry invariant if and only if k((px, ipy) = k(x, y) 
\/ip G Iso(H n ) in which case we say that k is doubly isometry-invariant. Working in the ball model 
and since every ip G Iso(H n ) is the composition of a translation with some U G 0(n), the additional 
requirement on the kernel 

k(x, y) =Y1 a i,j{Syx)w I (x) <g) w J {y) amounts to k(Ux, UO) = k(x, 0), that is, 

a/,j([/x)?7*u; / (x) <g> U*w J (0) = ^ a IyJ {x)w I {x) ® u> J (0), W 

Thus we are interested in describing those k(x,0) — which is a m-form at whose coefficients are 
m-forms in x — that are doubly invariant by all U G 0(n) in the sense above. Once the k(x,0) 
having this property are known, k(x,y) = k(S y x,0) defines the general doubly isometry invariant 
m-form. For m = the k(x,Q) are simply the radial functions a(|x|), and a{\S y x\) = a(\(p y x\) is 
the general doubly isometry invariant function. For m / their general expresion is not so simple. 
We find it more convenient to use the usual basis dx 1 so we look at k(x, 0) in the form 

k(x,0)= b IJ (x)dx I <g> dx J (0) (5) 

\I\=\J\=m 

and we must impose ^ bi t j{Ux)d(Ux) 1 ® d(Ux) J (0) = k(x,0) MU . For instance 



7 (x,0) = J2dx i ®dx i (0) 



is easily seen to be doubly 0(n)-invariant, and so is 



-L 7 A • • • A 7 = ^ d%1 ® dx\0) 



7m 

m! 

\I\=m 

(here we use the symbol A to denote as well the exterior product of double forms defined by 
(ai <g) Pi) A («2 <8> 02) = (cti A Q2) <8> A /?2))- Another doubly 0(w)-invariant 1-form is 

t(x,0) = ^Y^Xidx^ ® ^Y^Xidx^O^j . 

Lemma 2.1. The double forms 7 and r generate all doubly 0(n) -invariant k(x, 0). More precisely, 
their general expression in the ball model is 

k(x,0) = A 1 (\x\)j m + A 2 (\x\)T Aj m -! 0<m<n (6) 
k(x,0) = A(\x\)j m , m = 0, n 
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Proof. First we prove by induction the following statement S(n): if k(x,0) is a doubly invariant 
(p, g)-form cj 7 j dx 1 ®dx J (0) with constant coefficients, then k = if p ^ q, or k is diagonal 

|/|=P,|J|=9 

i.e. k(x, 0) = c dx 7 <g) dx 7 (0) = C7 P if p = q. Of course 5(1) is obvious; assuming S(n — 1), let 
\I\= P 

us break k(x,0) in four pieces, depending on whether = 1 or not: 

k= cijdx 1 ®dx J (Q)+ + 2 + 2 = f fc i + fc 2 + A; 3 + A; 4 . 

»i=ii=i ii=iji^l ii^i,ii=i u^ij'i^i 

We may write k\ = (dx 1 <g> <ix 1 (0)) A fei, fc 2 = (dx 1 (g> 1) A fe 2 , &3 = (1 <8> dx 1 (0)) A fe 3 , with fe 1; fe 2 , 
/c3, &4 double forms in the fix 2 , . . . , dx n , dx 2 (0), . . . , dx n (0) of bidegrees (p — 1, q — 1), (p — 1, q), 
(p,q — 1) and (p,q), respectively. Imposing that k is doubly invariant by U of the type 

/i o o -A 



Ui 



u 



Ui E 0(n - 1) 



(7) 



/ 



we see that fei, fe 2 , &3 and k^ are 0(n — l)-invariant. We apply the induction hypothesis: if p = q 
k2 = ka = 0, and fei, k^ are diagonal, i.e. 



k = ci V] dx 7 ® cte J (0) + c 2 V] dx 7 (8) dx / (0) 



(8) 



ii=l 



If we use now U E O(n) permuting the first two axes we see that c\ = c 2 and hence k is diagonal, 
establishing S(n) in case p = q. \i\p — q\ > 1 everything is 0. Finally if \p — q\ = 1, say p = q + 1, 
then fe 2 is diagonal and all others are zero 

k = c((ix 1 <g> 1) A Y dx ' J ® dx' J (0) 
\A=i 

where x' = (x 2 , . . . x n ). If we impose the invariance under the permutation of the first two axes as 
before, it is clear that k must be zero. 

Having proved that S(n) holds for all n, let now k(x,0) be as in (JSJ) doubly 0(n)-invariant. 
Clearly k(x, 0) is then determined by its values k(r , 0), where r = (r, 0, 0, . . . , 0). Fixed r, k(r , 0) 
may be regarded as a double (m, m)-form with constant coefficients, which is invariant by all 
U E 0(n) fixing r, that is, of type (J7J). We write now the decomposition of fe(r ,0) in terms of 
fei(r, 0), fe 2 (r, 0), &3(r, 0) and ^(r, 0) as before, and applying S(n) we get (JSJ) 

0) = ci(r) ^ dx 7 (7) ® dx 7 (0) + c 2 (r) ^ dx 7 (7) $ dx 7 (0) 

il=l il^l 
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(if m = n the last term is zero and the first is 7m ) which we write 

= (ci(r) - c 2 (r)) V] ® (fo 7 (0) +c 2 (r) ^ dx 1 (r) ® dx 1 ($) 

H=l \I\=m 
\I\=m 

= (ci(r)-c 2 (r))dx 1 (r)®dx 1 (Q)A ^ r) <g> dx T (0) + c 2 (r) ^ 7) tZx / (0) 

|/|=m-l |/|=m 

= (ci (r) -c 2 (r))r" 2 r(7, 0) 7m _i(7,0) +c 2 (r) 7m (7,0). 

Finally, fixed x we choose U such that Ux = r, r = |x|, and use the invariance of k, r, 7 to find 
©with A x (r) = c 2 (r), ^ 2 (r) = r~ 2 (ci(r) - c 2 (r)). □ 



To find the general expression of a doubly isometry invariant kernel k(x, y) we must translate 
k(x,0) to an arbitrary point: k(x,y) = k(S y x, S y y). We may use any isometry mapping y to 0, for 
instance we may use ip y given by © instead of S y . We introduce the basic forms a, (3, r and 7 

a = a(x,y) = <p*(x) dtp % y (x) , /3 = ^<P % v ip) &p\(v) = ~X^( X ) TZ1 12 

n 1 n 

r = a ® /?; 7(s, y) = ^ rf^(x) (8) = - _ ^ ^(x) ® d|/* = d x (3 . 

i=l i=l 

The lemma gives part (a) of the following theorem. Part (b) gives other equivalent general expres- 
sions, which are intrinsic, that is, independent of the model of M n at use. 

Theorem 2.2. (a) The general expression of an (rn,m)-form k(x,y) doubly isometry-invariant in 
M n , in the ball model, is 

k(x,y) = A 1 (\ip y x\)~/ m (x,y) + A 2 (\(p y x\) r(x,y) A7 m _i(x,y) , < m < n 

k(x,y) = A(\cpyx\)"f m (x,y) m = 0,n. 

(b) Another equivalent expression for < m < n is 

k(x,y)=B 1 (D)(d x d y D) m + B 2 (D)(d x D®d y D)A(d x d y D) m - 1 = 
= (d(D)d x d y D + C 2 (D)d x D ® d y D) m 

where D denotes an arbitrary function of the geodesic distance d(x,y). 

(c) All such k(x,y) are symmetric in x,y S EI n 

Proof. Part (a) has been already proved. For (b) note first that it is enough to consider one function 
of d: we choose D = r(x,y) 2 which in the ball model equals |(/7 y (x)| 2 . Then d x D = 2a, and using 
(HJ), © one finds 

d y D = 2(1 -D)J2 viWrAn* = ~ 2{ ~ l ~ D ^ • 
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This gives t = a® (3 



-l 1 

4 i-D^x-^ *v 1*2/ -l^, 



cL-D (8) cLD, and 



dadyD = +2d x D ® p - 2(1 - D)d x f3 = +4r - 2(1 - L>) 7 . 



Therefore (d x d y D) m and 2" 1 (1 — D) m 7m— l differ in a term containing r, and so (b) follows. 
Part (c) is a consequence of (b). □ 

We will need the expression of the generators r, 7 in terms of the invariant basis w l . We obtain 
these using formula (J2J) for r 2 (x,y) in the half-space model. First 



a 



d x r 2 



1 -r z 



n-l 



2 2(|x-y| 2 + 4x n2 / ri ) 



2 ^ x n (xj - yi)w l (x) + (2x n (x n - y n ) - |x - y| 2 )w n (x) 



i=i 



J r 2 



2(r 2_i) 2 (|x-y| 2 + 4x n y n ) 3 3 

In the following we write w 13 = w l (x) (g> w 3 (y). We have 

1 1 - r 2 

t = a® (3 



2 ^2vn(yj - Xj)m>{y) + (2y n (y n - x n ) - \x - y\ 2 )w n (y) . 



4 (\x - y\ 2 + 4x n y n ) 2 



1.1 



where the Pij(x,y) are certain homogeneous polynomials. As we know, everything can be written 
in terms of z = S,,x: for instance 



1 - r z 



4z n 



and say for i, j < n 



P; 



\x-y\ 2 + 4x n y n \z\ 2 + 2z n + I 



(\x - y\ 2 + Ax n y n ) (\x - y\ 2 + 4x n y n ) 2 (\z\ 2 + 2z n + l) 2 
Therefore we may write 



l-r 2 



(\z\ 2 +2z n + l) 2 



'■j 



For 7 = d x f3 we obtain a similar expression 



(9) 



n-l 



i,i=l 



2(xj - ^)(xj - yj) 
\x - y\ 2 + 4:X n y n 

/ n-l 



+ 



1 



2^ |xj - yi| 2 

i=l 

\x - y\ 2 + Ax n y r , 



n-l 



w n,n + £ 



2{xj - yj)(x n - y n ) 
^ \x-y\ 2 + Ax n y n 
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Again this can be written 

1 — r 2 \ — a „ , , „• „■ 1 — r 2 



l i3 



Notice that 

PiM = om q JiM = om 

(|z|2 + 2z n + l)2 U{ h (| z |2 + 2Zn + 1) 

and hence 

|r(x,y)| =0(l-r 2 ), | 7 (x, y)| = 0(1 - r 2 ) . (11) 

3 Riesz forms and Riesz form- potentials in HP 

3.1. Our next objective is now to find an explicit left-inverse L for A on D m (M n ). Since A is 
invariant by all isometries, L should be too. By what has been discussed in section 2, L should 
have a kernel k m (x, y), 

Lrj(x) = / r)(y) A * y k m (x,y) 



doubly invariant by all isometries. Alternatively, notice that if k is some kernel such that 

rj(x) = [ Arj(y) A * y k(x, y) , rj € V m (M n ) (12) 



(which formally exists because Ar] = 0, ij £ T> m (M n ) imply r] = 0) then its average over the unitary 
group 0(n) with respect the normalized left-invariant measure dfj,(U), 



ki(x,y)= / ka(Ux,Uy)dn(U) 

JO{n) 

still satisfies (|12() and it is doubly invariant by 0(n). If <p x is an isometry mapping x to 0, k2(x, y) = 
ki(ip x x,(p x y) is independent of <p x , satisfies 1)12(1 and is doubly invariant by all isometries. 

Anyway, we look for a doubly isometry-invariant kernel k m for which 1)12)1 holds, and then 
consider the operator L defined by k m as above. Taking for granted by now that this operator L 
is well defined on 2? m (HI n ) and maps P m (IH n ) into locally integrable m- forms, notice that 1)12)1 and 
the symmetry of k m together imply that L is a right-inverse too, that is, ALa = a for a S V m (Il n ) 
in the weak sense: 



(ALa, rj) = (La, At]} = / La(x) A * Arj(x) = 

J X 

a(y) A * y k m (x,y)^ A * Arj(x) 
<x(y) /\*y\ k m {x, y) A * Arj(x) \ = (a, rj) 



x (.Jy 
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We work in the ball model. By Theorem 2.2, k m (x,y) is of type 
km(x,y) = A(\(p x y\)^ m , m = 0,n 



k m (x, y) = Ai(\ip x y)^ m + A 2 (\ip x y\)r A 7 m _i , < m< n 

where 7 = ^d<p*(a:) ® d<p l x (y),T = a. ® (3 with a = ^ tp x (y)d(p x (x), j3 = ^(p l x (y)dip l x {y) (notice 

that we are exchanging x,y, using (c) in Theorem 2.2). Condition (|12|) implies A y k m (x,y) = 
in y 7^ x (while ALu; = w implies A x k m (x,y) = in x ^ y). In fact, (|T2*|) amounts to requiring 
A y k m (x,y) = 5 X in a sense to be described below. 



3.2. In a first step we look for conditions on the Ai, A 2 , so that A y k m (x,y) = in y ^ x. A 
lengthy computation will show that the general harmonic k m depends on four parameters. By the 
invariance of k m , we may assume x = 0, in which case, writing r = |y|, 

k m (x,y) = A(r)j m , m = 0,n 
k m (0,y) = Ai(r)7 m + A 2 (r)r A 7 m _i 

with 7 = dx r (0) (g> <iy*, r = a<8>/3,a = ^ y l dx l (0), (3 = rdr. Since k m (x, y) is again doubly 
invariant, it must have an analogous expression with m replaced by n — m. Indeed, it is easily 
checked that 

*x *v 1m = 7 \t(1 ~ r ) 1n—m 

[n — m)\ 

2\2m-n ( 2 In-m T A 7„_ m _i 



* y (rA^i) = (m-l)!(l-r 2 ) 2M ( 



(n — m)\ (n — m — 1)1 
whence * x * y k m (x, y) = ( n ™ m y_ (1 - r 2 ) 2m ~™7 n _ m for m = 0,n and for < m < n, 

( m — 1)1(1 — r 2\2m-n 

*x * y k m (0, y) = r ^ [(mil + r 2 A 2 )in-m - (n - m)i 2 r A 7n-m-i] • (13) 

(n — mj! 

Moreover, since * commutes with A, it is natural to require as well that * x * y k m = k n _ m , that is, 
we may assume from now on that < m < n/2. 
For m = 0, using ()3.3|) we find 

A(A(r)) = ^(1 - r 2 ) [-1(1 - r 2 )A" + ((3 - n)r + r~ x (l - n))A'] 

from which it follows that A'(r) = cq(1 — r 2 ) n ~ 2 r 1 ~ n and 



A(r) = ex - co Al - s 2 )"- V^" . 
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We start now computing A y k m (0, y) for < m < n/2 using that on m- forms A equals (— l) m+1 (* d* 
d + (—l) n d* d* ). The double form A y k m (x,y) is also doubly invariant, and therefore it must have 
the same expression as k m with A%, A 2 replaced by other functions B\, B 2 to be found. In the 
computations we will use besides (|13f) the equations 

d y a = >y, d y (r A 7 m -l) = —rdr A j m = —f3 A j m 



* y dr A7 m = (-l) r 



.(l_ r 2)2m+2-n r -l aA7n _ 



m—l 



(n — m — 1)! 

which are easily checked as well. First, d y k m (0, y) = (A\ — rA 2 )dr A j m , so by the equations above 



* x * y d y k m (0,y) =(-1) 



(1 - r 2 ) 2m+2 -™(^ - rA 2 )r~"a A 7n _ m _ 1 



(n — m — 1)|! 
def (-l) m m! 



(n — m — 1)! 



^3 « A 7, 



n— m— 1 j 



*xdy *« dyk m (0, y) 



{-l) m m\ 
(n — m — 1)! 



(^ 37n _ m + V A 7„_ m _i) 



* y ^ * a ^^(0, y) = (-l) m ( n m x ) * y * x (A 3 -/ n _ m + A' 3 r V A 7n- m -i) 



(-1) 



m(n— m— 1) 



m! 



2\n- 



(n — m — 1)! 



(1 - r 2 ) 



2n (n-m)!_ , , (n-m-1)! 



???-! 



7 m + A 3 r- 



???.! 



"7n 



,, _n (n — m — 1)! 



3 (m-1)! 

= ( _ ir (n-TO+l) (1 _ r2) „- 2m [ ((n _ + ^/^^ _ mA / r -l T A7m _ 1 ] . 

By the analogous computation, applying d y to (|T3*|) 

(m - 1)! 



*xd y * y k m (0,y) 



[(mA 1 +r 2 A 2 )(l 



2\n-2ml> 



(n — m)\ 

+ (n - m)rA 2 {\ - r 2 ) n ~ 2m ] dr A 7n _ m 

*y dy *y fc^O^) = (- 1) ("^H™"" 1 ) (1 - ^^+2-^-1 

= ~[(mA 1 + r 2 .4 2 )(l - r 2 )"- 2m ] ' + (n - m)r,4 2 (l - r 2 ) n " 2m ] a A 7m _! 
d|f (_l)(m+l)(n-m) A4aA7m _ i 

^ * y ^ * y fcm(0, y) = (-i)("-™)(™+i)(A' 4 rV A 7to-i + ^4 7m ) 
It follows finally that A = (-l) nm+1 (*d* d+(-l) n d* d*) on k m equals 

A y A: m (0, y) = Bu m + B 2 r A 7m _i 



(14) 



(15) 
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with 



Bx = -A 4 - (1 - r 2 )"- 2m ((n - m)A 3 + A' 3 r) 
B 2 = -A 4 r~ l +m(l - r 2 ) 7 



,.2 , ii 2ii/ j/ ^ I 



Therefore, A y k(0, y) = is equivalent to the system Si = 0, B 2 = 0. It easily follows from this 
that A3 satisfies the equation 

r(l - r 2 )A'l + [(n + 1) - r 2 (3n + 1 - Am)] A' 3 - 2(n - 2m)(n - m)rA 3 = . 

Replacing in the equation B\ = 0, A 4 by its expression in terms of A\ and A 2 , and then A 2 by its 
expression in terms of A\ and A3, we find that Ai satisfies the inhomogeneous equation 

r(l - r 2 )A'[ + [(n + 1) + (n - 1 - 4m)r 2 ] ^ + 2m(n - 2m)rA 1 = 

= 2rA 3 (l + r 2 )(l-r 2 ) n - 2m - 2 

The change of variables A\(r) = G(x), A%(r) = H(x), x = r 2 , transforms these into the hypergeo- 
metric equations 



x(l-x)H"(x) + 



71/ 3 

- + 1 - (-n + 1 - 2m)x 



77 

tf'(x) - (- - m)(n - m)ff = 



'71 71 1 77 

x(l - x)G"(x) + ^ + l - (2m + 1 - -)xj G'(x) - m(m - -)G 

= \{l + x)(l- x) n - 2m - 2 H(x) d ^ f(x) 



(16) 



This system is equivalent to A y k m (x,y) = in y / x, whence the general doubly-invariant A; m 
harmonic in y 7^ x depends on four parameters. Note that for m = § the homogeneous equations 
are the same and can be solved explicitely: the general solution is H = as~ 2 + 6 and 



G(x) = cx 2 



For m < t| a fundamental family for the first equation in (16) is given by 

m(x) = x 2F(-m, - - m, 1 - -,x),u 2 (x) = F(- - m, n - m, - + 1, x) 

The hypergeometric function in u\ is a polynomial in x of degree m with positive coefficients, 1 + x 
if m = 1. A fundamental family for the second is given by 

u 3 (x) = x 2F(m — n, m— -,1 - -,x) = x 2 (l - x ) n+1 F(— + 1 - m, 1 - m, 1 - — , x), 

n n 
u 4 (x) = F(m, m - -, 1 + -, x) 
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The hypergeometric function in u 3 is a polynomial of degree m—\ with positive coefficients (see 
|E] for all these facts). The wronskian w(x) for this second equation is, by Liouville's formula 

W(x) = W(x ) exp - f -2+ l -(^ + l --2) t dt = CmnX -f _ x n- 2m) 

Jxo t) 

It follows from this that the parametrization for G is given by 

G(x) = c(x)u%{x) + d(x)u4(x) (18) 
where c(x),d(x) satisfy, with H(x) = au\{x) + 6^2 (x), 

C ' (X) = x(l-x)W(x) = \ °™ H + X)XH1 ~ X) ~ 3MX) 

d'(x) = --^Mg- = -I c m 1 „if(x)(l + x)x* (1 - x)- 3 n 3 (x) . 
x(l — x)W (X) 2 

Once Ai(r) = G(r 2 ) and ^3(r) = H(r 2 ) are known, the kernel k m (x,y) is completely known, 
because by the definition of A 3 in (jUJ), A 2 (r) = -(1 - r 2 ) n " 2m_2 ^ 3 (r) + r~ l A[(r) = -(1 - 
x) n - 2m - 2 H{x) + 2G'{x). 

The choice a = 0, c(0) = (a = c = in the parametrization l|17[l for m = ^) gives all doubly 
invariant k m (x,y) which are globally harmonic, with no singularity, and they are therefore spanned 
by the forms corresponding to the choice G = U4 and to the choice a = 0, b = 1, c(0) = 0, d(0) = 0, 



G(x) = i / (l + t)(l-t)- 3 t n / 2 U2 (t) Ui (t)dt\u 3 (x) 



(1 + t)(l - t)~ 3 r /2 n 2 (t)u 3 (t) dt > u A {x) . 



As a particular case, note that for m = f ,7m is harmonic in H 2m , and it is the simplest example 
of a non-zero harmonic m-form in L 2 (H 2m ). 



3.3. Besides being harmonic in y ^ x, the singularity at y = x most be such that ()12j) holds. 
Again, we may assume x = 0; we check this property using second's Green identity, whose version 
for general forms we recall now. 

The operator 5 being the adjoint of d, one has for a smooth domain Q C B n and a, (3 smooth 
forms on O, with degreea = degree/3 — 1 

/ aA*P= daA*P- aA*5j3. 
Jon Jn Jn 

Given two m-forms rj, u>, applying this with a = 5rj, /3 = u, next with a = u>, (3 = dr\ and 
substracting one gets the first Green's identity for m-forms 

(5rj A*w-wA* drj) = / (Arj A*u — 5n A*5uj — dn A * duj) 
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Permuting uj,n and substr acting again gives the second Green's identity 



(5n A*iv — ujA*dn — 5ujA*n + nA* dw) = / (An A * cu — Aw A * n) . 
an Jn 

We apply this to Q = B(0, R) - B(0, e) < e < R < 1, n £ V m (W n ) and our k m (0, y) to get 



An A * y k m (0, y) = / (k m A * dn + 5 y k m A * n - 5n A * y k m - n A * dk m ) . (19) 
\y\>e J\y\=e 

In case m = 0, the terms in 5k m ,5r] are of course zero; to get a term in 77(0) on the right when 
e — > we need dk m of the order of e 1_n and k m of the order of e 2 ~ n in \y\ = e. That makes k m 
locally integrable too, and (|12|) is obtained leting e — > 0. This means that for m = k is unique 
and is given by the well-known Green's function 

A(r) = c n f (1 - s 2 )"- V-™ ds (20) 

J r 

for an appropriate choice of c n . In case m > 0, again we need \k m (0, y)\ = o(r 1 ~ n ) as r — > 0, so 
that the first and third terms on the right have limit as e — > 0; then k m is integrable in y and the 
integral on the left converges to J An A * k m . Using the expression for * dk m in ()14|) . we find 

/ r]A*dyk m = — — — A 3 (e)* x n A a A 7 n -m-i ■ 

y|v|=e (n-m-1)! 

By Stoke's theorem, and since a = 0(r), the last integral equals 

(-l) m / A 7 n- m + O(e) . 

If ^3(5) = aoe _n + . . . , we see that 

lim / 77 A * dy/c m = c n (n — m) ml ao^(0) . 
£ J\y\=e 

Using (|15() for <5/c m = (— l) n ( m + 1 )+ 1 * d* , and proceeding in the same way, 



5 y k m A *n = -A 4 (e) / aA7 m _iA*i) 
'\y\=e 



= -A 4 (e) / ( 7m A* n + 0(e)). 
J\y\<e 

But by the equation B x = 0, A 4 (e) = -(1 - e 2 ) n ~ 2m ((n - m)A 3 (e) + eA' 3 (e)) = a me- n + 0(e 1 ~ n ), 
and hence the limit of the above expression is — c n m\cLQmn(0). Altogether, we conclude that if 
^3(5) = ao£~ n + o(e 1 ~ n ) and k m (0,y) = o(r 1_ "), one has 

An A * y k m (0, y) = -c n n ml a n(0) 
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so (|12|) will hold for an appropriate choice of a^. Taking into account the definition of A3 in (|14|) 
and that \k m \ ~ \A\\ + r 2 !^!; we see from (|17|l that if m = § this is accomplished by the choice 
c = 0, a = a ; then G(x) ~ logx, A 1 (r) ~ logr, A' x (r) = 0(l/r), A 2 (r) = 0(r 2 ) if n = 2; if n > 2 
,^4i(r) ~ r 2 ~ n anc l ^4 2 = 0(r~ n ). For < m < ^, in terms of the functions H, G introduced before, 
this translates to H(x) ~ aox" , G(x) ~ x 1 " 2". Now look at the general expression of H, G in (|18|) . 
The condition -ff(x) ~ co x ~^ fixes a = ao; then near x = c'(x) is bounded and d'(x) behaves like 
x~2. Since U4(x) is bounded, the term d{x)ui{x) behaves like x 1_ 2. So, we must normalize c(x) 
by c(0) = 0, so that c(x) = 0(x) and the other term c{x)u^{x) will behave like x 1_ 2 . 

In conclusion, all this discussion shows that the doubly invariant kernels k m (x,y) satisfying 
(fT2|) constitute a two parameter family described by H = a$ui(x) + bu 2 (x), c(0) = 0. The two 
parameters are b and the constant of integration for d(x) in 1)18(1 . Equivalently, they are obtained by 
adding to the form corresponding to H = aoui(x), c(0) = and say = the general globally 
smooth one described before. 



3.4. In order to produce the best estimates, in a sense we need to choose the best of the kernels 
k m . Naturally enough, we choose the k m having the best behaviour at infinity, x = 1, that is, 
so that G, H have the best decrease in size as x — > 1. In case m = ^, where we already have 
the normalization c = 0, a = ao, the choice b = —a gives the best growth H(x) = 0(1 — x) and 
G{x) = 0(log(l -x)). 

The hypergeometric function 113 behaves like (1 — x) n+1_2m near x = 1 while u^{x) = F(m,m — 
§, 1 + §, x) is bounded because 1 + ^— m — (m — ^) = 1 + n — 2m > 0. Similarly, u\ is bounded 
near x = 1; for U2(x) = — m, n — m, ^ + 1, x) we have ^ + 1 — — m) — (n — m) = 2m + 1 — n 
and hence it behaves like (1 — x ) 2m+l ~ n if 2m < n — 1 and like log(l — x) if 2m = n — 1. We use 
equations (|TH)) 

C{x) = C myn [ H{t)(l + t)t"/ 2 (l - t)~ 3 Ui (t)dt 

Jo 

d{x) = -c m>n [ H(t)(l + t)t n ' 2 {l - t)- 3 u 3 {t)dt + d . 

If b ^ 0, then = a «i(t) + 6n 2 (t) behaves like (1 - t) 2m+1 ~ n if 2m < n - 1 and like log(l - £) 
if 2m = n — 1, resulting in c(x) = 0(1 — x) 2m_n_1 , <i(x) = 0(log(l — x)) if 2m < n — 1 and 
c(x) = 0((1 — x)~ 2 log(l — x)), d(x) = 0((1 — x)" 1 log(l — x)) if 2m = n — 1. So if b / one has 
G{x) = 0(log(l-x)) if 2m < n-1 and G(x) = O^l-x)- 1 log(l-x)) if 2m = n-1. If b = 0, then 
ii" is bounded, giving c(x) = 0((1 — x)~ 2 ) and d(x) = 0(1) for 2m < n — 1, d(x) = 0(log(l — x)) 
for 2m = n — 1. In case 2m < n — 1, however, we can choose the constant do so that d(l) = 0, 
and then d(x) = 0(1 — x) n ~ 2m ~ x . This choice gives G(x) = 0(1 — x) n_2m_1 for 2m < n — 1. For 
2m = n — 1, no choice of do can improve the bound G(x) = 0(log(l — x)). 

It remains to estimate the growth of A 2 (r) near r = 1. Recall that the definition (|14() of A3 
translates to A 2 (r) = 2G'(x) - (1 - x)"~ 2m ~ 2 iJ(x). Both terms growth like (1 - x) n ~ 2m " 2 , but a 
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cancellation occurs. The functions 111,113 are C°° at 1 and have developments 

u 3 (x) = A(l - x ) n+1 ~ 2m + 0(1 - x ) n+2 - 2m 

u' 3 (x) = -A{n + 1 - 2m) (1 - x) n " 2m + 0(1 - x) n+1 - 2m 

H(x) = a ui(x) = B + 0(1 - x) . 

In Ui(x) = F(m, m — ^, 1 + ^, x), 1 + ^ — m — (m — ^) = n + 1 — 2m > 2, whence U4 has a finite 
derivative at 1 and a development 

u 4 (x) = C + D(l - x) + 0{1 - x) 1+£ Ve<l, u 4 (x) = 0(1) . 

Then W(x) = u 3 u 4 — u 3 u 4 = CA(2m — n — 1)(1 — x) n_2m + . . . , and so the constant c mn in (fTH|) 
is CA(2m - n - 1). Then from ((IHJ) 

c ' (x ) = ^a-^r 3 - + Oil - xy 2 

v 1 A(2m-n-\) y 1 

d'(x) = w + 0(1 - x)- 2 - 1 

C (2m — n — 1) 

which gives 

c(x) = -(1 - x)- 2 + 0(1 - x)" 1 

v ; 2 2(2m - n - 1) 

d(x) = 0(1 - x)™" 2 " 1 " 1 , 2m < n - 1 , 0(log(l - x)), 2m = n - 1 . 

But O' = c(x)ii' 3 (x) + <i(x)u 4 (x); the second term d(x)u' 4 (x) satisfies the required bound while the 
first c(x)ii' 3 (x) has a development 

c(xK 3 (x) = 7^ ~aTo ~ U A ^ n + 1 - 2m)(l - x)- 2m - 2 + 0(1 - xf- 2 "- 1 

2 A(2m — n — 1) 

= ^(1 _ x )«- 2m - 2 + 0(1 - x) n - 2m - 1 . 

As (l-x) n - 2m ~ 2 i?(x) = J B(l-x) n - 2m - 2 + 0(l-x) n ~ 2m " 1 the bound for A 2 follows for 2m < n-1. 
However, for m = this no longer holds. Indeed, from (|T7|) . where c = 0, a = ao, 6 — a) 

2G'(x) = x-f" 1 r s"/ 2 (l + s)(l - s)- 3 a( S -% - 1) 

has development 

2G'(x) = na(l - x)" 1 + 0(log(l - x)) 

while 

(1 - x)" 2 a(x"? - 1) = |a(l - x)" 1 + . . . 
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We point out that all this can be obtained, in loose terms, working directly with the hypergeo- 
metric equations relating G, H 

x(l - x)G"(x) + |J + 1 - (2m + 1 - |)x] G'(x) - m{m - ^)G = i(l + x)(l - x) n ~ 2m ~ 2 H{x) 

and using asymptotic developments. If H(x) = ho + h\(l — x) + . . . and G(x) = gj(l — x)- 7 + . . . , 
identifying the lower order terms in both sides gives, 

9j j(j - 1 - n + 2m)(l - x)^ 1 = fto(l - x)"" 2 " 1 " 2 . 

When H = 0, one must have either j = (corresponding to U4) or j = n-2m+l (corresponding to 
U3). For the inhomogeneous equation, if j 7^ 0, j 7^ n + 1 + 2m (that is G contains no contribution 
from U3, U4) one finds j = n — 2m — 1 if 2m < n — 1 and gjj = — 4p. Then 2G'{x) = h${\ — x ~) n ~ 2m ~ 2 + 
. . . , (1 — x) n ~ 2m ~ 2 H(x) = ho(l — x) n ~ 2m ~ 2 , showing cancellation. An analogous argument works 
if 2m = n — 1, but not for 2m = n. 

We summarize the results in this and the previous subsections: 

Theorem 3.1. For \n — 2m| > 1, there is a unique doubly invariant kernel 

km(x,y) = Ai(\ip x y)j m + A 2 (\<p x y\)T A 7 m -i,m 7^ 0; k m (x,y) = A(\(p x y\)j m , m = 0,n 

/or which (|12l) holds, and satisfying moreover 

\Ai{r)\ = 0{l - r 2 )\ n ~ 2m \~ l as r -» 1 . 

For m = , i/iere is a one-parameter family of such kernels satisfying 

\A i (r)\ = 0{\og{l-r 2 )). 

For m = 7j , there is a one-parameter family of such kernels satisfying 

\A t (r)\ = 0(1 - r 2 )- 1 . 

In all cases A\(r) ~ r 2 ~ n , A 2 (r) ~ r~ n as r — ► 0. 

For |n — 2m I > 1, we call k m (x,y) the Riesz kernel for m-forms in BP, and 

Lr)(x) = rj(y) A * y k m (x,y) 



the Riesz potential of 77, whenever this is defined. From IjllJ) we see that 

Mz,y)|=0(l-r 2 ) n - m - 1 . (21) 

With the notations used before, the function A^(r) = H[r 2 ) is bounded with bounded derivatives 
near r = 1. Then (|14|) and symmetry imply 

KMx,y)|, |d y A; m (x,y)| =0(l-r 2 )"- m - 1 (22) 
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too. The growth of ^3 also implies A4 = 0(1 — r 2 ) n 2m because B\ = 0, and then gives as 
well 

\S x k m (x,y)\, \5 y k m (x,y)\ = 0(1 - r 2 )™-™- 1 . (23) 

By construction, one has LAw = n for r] E T> m (M. n ). We will need the following generalization 
of this fact. 

Proposition 3.2. If r/ is a smooth form in M n such that 

\v(y)\,\^r ] (y)\=o(l-\y\ 2 r, y G B" 

then LAn = r]. 

Proof. In (|19|) we would get an extra term 

(k m A * dr] + 5k m A*n — 5n A * k m — n A * dk m ) . 



'\y\=R 

Estimates dHJ), (jUJ) and (J23J) imply that with x fixed and \y\ = R / 1 

Inserting |ry(y)|, |V?7(y)| = o(l — |y| 2 ) m we see that this extra term vanishes as R f 1. □ 

4 Proof of the main theorem 

4.1. Once the Riesz form k m (x, y) has been found, our aim is now to prove that the corresponding 
convolution 

L m rj{x)= \ r/(y) A* y k m (x,y) 
Jh» 

satisfies 

||L m 77||p )S+2 < c||7?|| P)S (24) 

for m 7^ and p in the range pi(m) = w " ] _ 1 , m < p < = P2{m), and for a compactly 

supported m-form n (recall that we are assuming without loss of generality that m < ^). Since 
these are dense in the Sobolev spaces and we already know that AL m rj = L m Arj = n, this will 
prove the theorem for m 7^ ^w^, §• The case m = will be commented later. 

We work in the translation invariant basis w 1 . Taking into account formulas Q and (10) for 
7, r, the Riesz- form is written in the model 

k m (x,y)= ^2 a I)J (Syx)w I {x)®w J (y) 

\I\=\J\=m 

where each coefficient aij has an expression, with z = S y x 

f \ lT , / \ PlA z ) 2 l + M 2 -22; n \x-y\ 2 



(\z\ 2 + 2z n + l) 2m ' l + |z| 2 + 2z n |x - y\ 2 + 4x n2/ri 
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Here pij(z) is a certain polynomial in zi,...,z n , ^ij is C°° in (0, 1) with ^/^(r) ~ c^r 2 n as 
r \ 0, */,j(r) = 0(1 - r 2 )™-" 1 " 1 as r / 1. The term qj,j(z) = Pi,j(z) / '{\z\ 2 + 2z n + l) 2m is 
bounded. 

If f? = i the coefficient (Lrj)i(x) of -L77 in the basis w 1 is a finite linear combination 

of hyperbolic convolutions 

(Lr]) I (x) = V / ^i,j{r)qi t j{z)r]j{y) dfi(y) . 
j J ^ n 

By ellipticity of A, Lrj is a smooth form. Moreover, since 77 has compact support, we see from (j2J) 
and (EHJ) that, in the ball model, 

\L V (x)\, \d(L V )(x)\, \5(Lrj)(x)\=0(l-\x\ 2 r- m - 1 

which amounts to 

\(L V )j(x)\, iXiiLr,)^)] = 0(1 - \x\ 2 r- m ^ . (25) 
We claim that for second-order derivatives we have too 

{XjXiiLri)^ = 0(1 - \x\ 2 ) n - m - 1 , i.e. \V^{Lrj)(x)\ = 0(1 - M 2 )™-" 1 - 1 . (26) 

Notice that since we already know that AL77 = 77, from the expression of A in the basis w 1 given 
in (3) it follows that it is enough to show that for j < n. We will see below (equation (30) and 
invariance of the Xj) that each of the functions a(z) = VP i 1 j(r)qi t j(z) satisfies 

\X 3 X ia {z)\ = 0(1 - r 2 ) n - rn ~ l 

from which (|26|) follows as before. In fact, the discussion that follows will show that \\7^ Lrj(x)\ = 
0(1 - \x\ 2 ) n ' m - 1 VJfe. 

We continue the proof of (|24jl. We claim first that it is enough to prove (|24jl for s = 0. For 
a smooth form 77 = ^rijw 1 let JQ77 denote here the m-form JQ77 = ^Xirijw 1 . It is clear from 
formulas (3) and the commutation properties 

[Xi,Xj] = , i,j <n, [X 1 ,X i ]=X i , i<n 

that for each i there is an operator Pj of order two in the Xi , . . . , X n such that 

XiAr} - A(Xiri) = Pi{X) V . 

Applying this to L77, which is smooth by the ellipticity of A, we get (Xj — AXiL)r] = Pi(X)Lrj. 
But XiLrj satisfies, by and (gHJ) 

|XLt7(x)|, \d{XiLn){x)\, \5{X i Li 1 ){x)\=0{l-\x\ 2 ) n - m - 1 

and hence by Proposition 3.2, LA = Id on it. We conclude that for all 77 S 2? m (BI n ) 

{LX i -X i L)r} = LP i {X)U l . 
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Assume that (|24jl has been proved up to s, so that by density it holds for a £ H^^W 1 ) too, and 
let 7 be a multiindex of length | —y | < s. For i = 1, . . . , n and r] £ D m (W n ), 

X^XiLr] = X^LXtT] - X"<LPi(X)Lrj 

so using twice the induction hypothesis 

\\X^XiLr,\\ p < const (\\X ir] \\ P)S + \\Pi(X)L V \\ PjS ) < 

< const (||r/|| PiS+ i + \ \r]\\ PtS ) 

proving (|24j) for s + 1. Proving (|24j) for s > means proving 

||(L77)/|| P , ||Xi(L77)/|| p , ||XjXi(L7?)/|| p < const | | p . 



As before, using that we already know that ALr/ = rj we see that for the second-order derivatives 
we may assume j < n. In the following we delete the indexes /, J and denote by a{z) = ip(r)Q(z) 
a convolution kernel with ip, Q as above, and proceed to prove that the convolution 



(C a a)(z) = a(S y x)a(y) d/J,(y) 
Jw n 

satisfies 

||C a a|| p , \\Xi(C a a)\\ p , \\XjXiC a (a)\\ p < const \\a\\ p , pi < p < p 2 (27) 

where in the last case we may assume that j < n. The fields Xi are invariant, and therefore 
XiC a a, XjXiC a a are obtained, respectively, by convolution with Z^a, ZjZiCt (by (jU). Recall that 

(\ n— m— 1 
i+\$+2z n j as r / 1 and V(r) ~ r 2 ~ n as r \ 0. 

In order to estimate Z^a, ZiZja, we collect first some auxiliary estimates. We claim that 



\ZiQ\ < const, \ZiZjQ\ < const 

\Z%r\ < const(l — r 2 ) , \ZiZjr\ < const r _1 (l — r 2 ) . 
The first two are routinely checked, for instance, when differentiating the denominator in Q, 



(28) 



1 



(1 + \z\ 2 + 2z n ] 



2m 



4mz n zi 



(l + |z| 2 + 2z n ) 2m + 1 



const 

~ (l + |z| 2 + 2z n ) 2m 1 ; 



so that the term pi y j(z)Zi [(1 + \z\ 2 + 2z n ) 2m ] will still be bounded. All other terms can be 
treated similarly. Differentiating 1 — r 2 = 1 , %" „ — we get 

_l-r 2 Zi z n l-r 2 ll + |z| 2 -2zl 

6i.T — t. ; — tt; ; r , ^JnX 



2 r(l + \z\ 2 + 2z n ) ' 2 rl + \z\ 2 + 2z r . 



77 - 1 ~ T<1 f 5 jl*k 1 + z i z 3 z n ( ■ 

3 ir ~ 2r \\ + \z\ 2 + 2z n 2r 2 (1 + |z| 2 + 2z n ) W ' '" J """ 



j ™ r " r 1 (l + |z|2 + 2 z n )2) + 4r 2 (1 + \z\ 2 + 2z n f ' ' 7 ' " ' 
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These imply (|2*H|) because 

1*3*1,1 + \z\ 2 - 2z 2 n < (1 + \z\ 2 - 2z n ) l / 2 {l + \z\ + 2z n ) 1 ' 2 = r(l + \z\ 2 + 2z n ) . 

Now 

Z ia {z) = i>'{r)Z i rQ{z)+^{r){Z i Q){z) 

ZjZ ia (z) = i)"{r){Zir){Z jr )Q{z) + tf(r){ZjZir)Q + ^'(r^r^Q 
+ i;'{r)Z jr ZiQ + ^(rjZ^Q . 

The estimates (|2*H)l imply 

|a(z)|,|Z i a(z)|,|2' i Z i o(z)| = 0(l-r 2 ) n - m - 1 as r/l 
\a(z)\ = 0(r 2 - n ) , \Z ia (z)\ = 0^) , \Z 3 Z t a{z)\ = 0(r~ n ) 



(29) 



(30) 



We will call a convolution kernel b(z) m- admissible if \b(z)\ = 0(r 1 ~ n ) as r \ and, moreover, 
|6(z)| = 0(1 — t- 2 )™-™- 1 as r y 1. We will prove later (Theorem 4.2.) that hyperbolic convolution 
with m-admissible kernels defines a bounded operator in L p (EI n ) for the range (pi(m) , p2(m)) 
specified in the statement of the main result. From the estimates (|29|) we see that a and Zja are 
m-admissible kernels, and so (|27[) will be proved for them. As \ZjZ^a{z)\ = 0(r~ n ) has the critical 
non-integrable singularity at r = 0, ZjZia(z) is not an m-admissiblekernel. Notice however from 
(|29|). (|3U|) that the last three terms ip' (r)ZirZjQ, (r)ZjrZiQ, ip(r)ZjZiQ are indeed m-admissible. 
Moreover, the estimate \ZiQ\ < const implies that Q is Lipschitz with respect the hyperbolic metric, 
in particular 

Q(z) = Q(e) + 0(log ^±1) = Q(e) + O(r) 
1 — r 

for small r. This means that replacing Q by Q — Q(e) in the first two terms leads to an m-admissible 
kernel again. All this leaves us with the kernel 



ip" (r)ZirZjr + ip (r)ZjZir , j < n . 

If ip(r) = cor 2 ~ n + . . . , write <p(r) = cor 2_n (l — r 2 ^ n ~ m - 1 ^ then the above differs from 

(f>" (r)ZirZjr + 4>'(r)ZjZir 

in an m-admissible kernel. By the same reason, we may replace (j>"(r), 4>'{r) respectively by 
(r 2-n )"(l — r 2 ) ri ~ m ^ 1 5 (r 2 ~ n )'(l — t- 2 )™-™-^ that is to say we must deal with the convolution 
kernel 

(1 _ r 2 ) n " m - 1 Z j Z i (r 2 " ri ) . (31) 

We introduce a class of singular hyperbolic convolution kernels to deal with the later. For this 
purpose it is more convenient to work in the ball model, so now b is defined in B™ and r = \z\. We 
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replace the integrable singularity r 1 n by a typical Calderon-Zygmund singularity (see e.g. [3]). 
Thus, we will call b a m- Calderon-Zygmund singular kernel if it has the form 

b(z) =n(w)r- n (l-r 2 ) n - m - 1 , z = rw, w G 

where f2 is say a Lipschitz function on S" -1 satisfying the cancellation condition 

Sl(w) da{w) = O . (32) 



In Theorem 4.2. below we prove that m-Calderon-Zygmund singular kernels define bounded oper- 
ators in the same range of p. With the following proposition, applied to <f>% [z] = |z| 2_ra this will 
end the proof of the main result. The proposition is the analogue of the well-known statement that 
for (j) smooth and homogeneous of degree 1 — n in R n , ^ defines a Calderon-Zygmund kernel (it is 
homogeneous of degree — n and the cancellation condition (32) is automatically satisfied, because 
/n<M<m & dV ^ = (Jl«|=m " h\=n) ^ A • • • A dx i A • • • A dx n = 0). 

Proposition 4.1. If(f>i,<j>2 are homogeneous functions of degree 1 — n, 2 — n respectively, the kernels 
(1 — r 2 ) n ~ rn ~ l Zi(f>\, (1 — r 2 ) n_m_1 ZjZj</>2 are sum of (m — l)-admissible and (m — 1)- Calderon- 
Zygmund singular kernels. 

Proof. We replace the Zj by Yj = (1 — r 2 )-^-; we have 



Yifa) = (1 - r 2 )^i y f 0a = (1 - r 2 )^ = (1 - r 2 )0{r 



so in all cases we get an extra factor (1 — r 2 ). Besides q z .q 2 z . are, as noted before, homogeneous 
of degree — n, and satisfy the cancellation condition (|3*2"|) . □ 



4.2. It remains to prove 

Theorem 4.2. Both m-admissible and m- Calderon-Zygmund kernels define by hyperbolic convo- 
lution bounded operators in L p (M. n ) for , < p < 2—^, < m < 

We will make use of the following well-known Schur's lemma for boundedness in IP of an integral 
operator with positive kernel. 

Lemma 4.3. // K (x, y) is a positive kernel in a measure space X and 1 < p < oo, the operator 
Kf(x) = f x K(x,y)f(y)dfi(y) is bounded in L p ([j) if and only if there exists h>0 such that 

[ K{x,y)h{yf d^y) = 0{h{xY) , xeX (33) 
Jx 

[ K{x,y)h(x) p d^x) = 0(h(y) p ), y£Y. (34) 
Jx 
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Here q is the conjugate exponent of p, - + ~ = 1. If h can be taken = 1 that is 

sup / K(x, y) dfj,(y) , sup / K (x, y) dfi(x) < +co 
x Jx y Jx 

then K is bounded in L p {jj) for all p, 1 < p < oo. 

Proof. Let us prove Theorem 4.2. If b is m-admissible, b = 61 + 62 with 61(2) = 0(r 1_n ) for r < 1/2, 
61 (z) = for r > 1/2, and 62(21) = 0(1 — r 2 ) n - m_1 for all r. We apply to b\ the second criteria in 
Lemma l4.HI working in the ball model (recall that 15^X1 = \(f y x\ is symmetric in x,y) 



/ bi(Syx) d/j,(x) , / bi(Syx) dfi(y) < c / l<Sy x | 1 n d/J,(x) 

JX JX J\S y x\<\/2 

f 1 ll-n^ / \ /" 1/2 
c / z cmz = const / ^— < +00. 

J\z\<l/2 ^ ' Jo (1 " r 2 )" 



We apply to (1 — r 2 ) n m 1 the criteria of the first part on Lemma 14. 3( working this time for 
convenience in the half-space model, where the kernel is written 

n—m—l / n \ n—mi 

K(x,y) = (l-r 2 ) r 



2\n— m— 1 / ^ z n \ ( 4x n y n 



1 + \z\ 2 + 2z n J \\x - y\ 2 + 4:X n y n 

We test h(y) = y" in for an exponent a to be chosen, so we need 



r —k—l+aq j 

/ Vn a V _ q faq+k+l-n 

Jy n >o (\x ~ y\ 2 + Ax n y n ) n - k - 1 ~ 



We write \x — y\ 2 +4:X n y n = \x' — y'\ 2 + (x n +y n ) 2 , where x' = (x\, . . . , x n -±) £l" 1 and analogously 
for y' , and integrate first in y 1 . One has for 2m < n — 1 



/ 



dy' f°° s n ~' 2 



(\x' - y'\ 2 + (x n + yn) 2 ) 71 -™- 1 Jo (s 2 + (x n + yn) 2 )"- 1 

= 0{(x n + y n ) 2m+l - n ) 

and so the above becomes 

,oo y « g -m- l dyn 



lO \X n 1 Vn) 

By homogeneity (y n = x n t) this reduces to 

poo ^a.q—m—1 



0(x 



aq+m+l— n\ 
n ) 







[l + t 



in— 1— 2m 



O(l) 



which holds whenever m < aq < n — 1 — m. By symmetry, for (|34|) we need as well m < ap < 
n — 1 — m. Therefore, a choice of a is possible whenever mmax 1 ^ ' • 1 " 1 ^ 1 1 
and this gives the range n ™i} m < p < ■ 
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Consider now a m-C alder on- Zygmund kernel b(z) = Q(w)r n (l—r 2 ) n m , Since \S y x\ = \(f y x\, 
we may replace z = S y x hy z = ip x y. Using Q this is given by 

(x — y)(l — \x\ 2 ) + x\x - y\ 2 
Z ~ A 

where we use the notation A = (1 — |x| 2 )(l — |y| 2 ) + \x — y\ 2 ; note that 

(i-M 2 ),(i-M 2 )<^ 1/2 

Also recall that r = \z\ and Ar 2 = \x — y\ 2 . Hence we can write 



z x — y x — y f I — \x\ 2 



1 + x ■ r . 



r \x-y\ \x-y\\ y/A 
But 

1- \x\ 2 _ i _ (1 - \x\ 2 ) 2 -A (1- \x\ 2 )Q(\x - y\) + Q(\x - y\ 2 ) 

VA ~ VA({1 - \x\ 2 ) + y/A) ~ A 

is 0(r). Therefore, modulo an m-admissible kernel, we may replace Q(w) by $7 ^ \x^y\ ) ■ This leaves 
us with the kernel 

K = (l- r 2 ) n ~ m - l VL ( ^—^) r -n = n_ ^n-m-l^ _ n ( IJZJL) A % (x, y) . 

\\x-y\J \\x-y\J 

Fix p, 1 < p < oo. Write 

A2( x ,y) = (1 - \xf)p(l- \y\ 2 )i +0[\x-y\A~\ 

Since \x — y\ l ~ n A~^~ = r 1 ~ n , the kernel K differs from 

(1 _ r 2 ) n - m - l \x - y\- n n (1 - |x| 2 )?(l - \y\ 2 f« 



\x - y\ 

in a m-admissible kernel, so we keep this one. We write it as the sum of 



\x-y\~ n n {j^fj (1-N ) 7 =Ki(x,y) 

and another K2(x,y) which we estimate by 

\K 2 (x,y)\ = O (r 2 \x - y\~ n (l - |x| 2 )?(l - |y| 2 )?) = O (r 2 ^(l - |x| 2 )?(l - \y\ 2 )^A~^ . 

Write K^i for the (euclidean) C alder on- Zygmund convolution operator with kernel 
\ x ~ y\~ n ^ (jf^ff)' wn i c h as h is well-known, satisfies an L p (dV )-estimate. Notice that 

K x f{x) = (1 - \x\ 2 )^K n (/(l - \y\ 2 r?) 
and therefore, using the L p -boundedness of 

\K 1 f(x)\Pd fi (x)= [ K n (f(l - MT*) P dV(x) < [ \f(x)\'dn(y). 
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For K2, we can ignore the integrable simgularity r 2 n and arguing as we just did with K\, we need 
to show that the integral operator 

K ^m= [ n 1 rri ^f(.y)dv( y ) 

J\ y \<i (1 - \x\ + \x- y\) n 

satisfies L p (dy)-estimates for all p, 1 < p < 00. To see this, just check that the criteria in Lemma f4.3l 
holds with h(x) = (1 — |x| 2 ) ~pq 

□ 

Notice that in case m = a m-Calderon-Zygmund kernel defines a bounded operator in all 
L p (M n ), 1 < p < 00: this is the right analogue of the euclidian kernels, because (1 - r 2 ) n ~ l is the 
typical growth at infinity of a weak L 1 (d/i) function in HP. 

4.3. Finally we make some comments, with no proofs, on the critical case m = rj ^- in the 
main Theorem. In this case the m-admissible and m-Calderon-Zygmund operators appearing in 
XjXiC a u, etc. have (1 — r 2 )^ log instead of (1 — r 2 Y~ m ~ l = (1 — r 2 ) 2 ^ as a factor. One 
can then prove that for (3 > and 2 < p < 2 + fa—i) ' 

1 1 V?l La < const / |77P(1 - ly] 2 )^ dpi(y) . 

The L p -estimates do not hold in this case for any p, because they do not hold for p = 2 and A is 
self-adjoint. 
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